AgRISTARS 


'1041 

SM-T3-04436 
NAG 5-31 


A Joint Program for 
Agriculture and 
Resources Inventory 
Surveys Through 
Aerospace 

_ .... . Remote Sensing 

Soil Moisture 

August 1983 


Technical Report 

A MATHEMATICAL CHARACTERIZATION OF 
VEGETATION EFFECT ON MICROWAVE 
REMOTE SENSING FROM THE EARTH 


Y. Choe 
L. Tsang 





1. Report No. 

SM-T3-04436 


2. Government Accession No. 


4 Title and Subtitle 

A Mathematical Characterization of Vegetation Effect on 
Microwave Remote Sensing from the Earth 


3. Recipient's Catalog No. 


5. Report Date 

August 1983 


6. Performing Organisation Code 


7. Author(s) 

Y. Choe 
L. Tsang 


9. Performing Organization Name and Address 

Remote Sensing Center 
Texas A&M University 
College Station, Texas 77843 


12. Sponsoring Agency Name and Address 

National Aeronautics & Space Administration 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 


.16. Abstract 


8. Performing Organization Report No. 

RSC-139 


10. Work Unit No. 


11. Contract or Grant No. 

NAG 5-31 


13. Type of Report and Period Covered 

Quarterly Report 


14. Sponsoring Agency Code 



In passive microwave remote sensing of the earth, a theoretical model that 
utilizes the radiative transfer equations has been developed to account for the 
volume scattering effects of the vegetation canopy. 

Vegetation canopies such as alfalfa, sorghum, and corn are simulated by a 
layer of ellipsoidal scatterers and cylindrical structures. The ellipsoidal scat- 
tered represent the leaves of vegetation and are randomly positioned and orient- 
ed. The orientation of ellipsoids is characterized by a probability density func- 
tion of Eulerian angles of rotation. The cylindrical structures represent the 
stalks of vegetation and their radii are assumed to be much smaller than their 
lengths. The underlying soil is represented by a half-space medium with a homoge- 
neous permittivity and uniform temperature profile. The radiative transfer equa- 
tions are solved by a numerical method using a Gaussian quadrature formula to 
compute both the vertical and horizontal polarized brightness temperature as a 
function of observation angle. 

The theory was applied to the interpretation of experimental data obtained 
from sorghum covered fields near College Station, Texas. 


17. Key Words (Suggested by Author(s) I 


18. Distribution Statement 


Unlimited 


19. Security Qassif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21. No. of Pages 
107 


22. Price' 


'For sale by the National Technical Information Service. Springfield. Virginia 22161 


JSC Form 1424 Irtev Nov 75) 


NASA JSI 

















SM-T3-Q4436 


NAG 5-31 


A MATHEMATICAL CHARACTERIZATION OF VEGETATION EFFECT ON 
MICROWAVE REMOTE SENSING FROM THE EARTH 


by 

Y. Choe 
L. Tsang 


This report describes activity carried out in support of the Soil 
Moisture activites of the Project. 


Remote Sensing Center 
Texas A&M University 
College Station, Texas 77843 


August 1983 


TABLE OF CONTENTS 


Page 

INTRODUCTION 1 

DESCRIPTION OF MEASUREMENT AND VEGETATION 8 

Background 8 

Data Acquisiton 8 

Description of Vegetation 11 

RADIATIVE TRANSFER THEORY 14 

Radiative Transfer Equations 15 

Boundary Conditions 21 

Cylindrical Scatterer Model 23 

Formulation 23 

The Solution for Cylindrical Model 26 

The Results 27 

Ellipsoidal Scatterer Model 38 

Formulation 38 

The Solution and the Results 44 

Composite Form of Cylindrical and Ellipsoidal 

Scatterer 58 

NUMERICAL APPROACH 60 

DATA MATCHING 67 

CONCLUSION 82 

REFERENCES 84 

APPENDIX A 89 

APPENDIX B 92 

APPENDIX C 100 

APPENDIX D 104 


LIST OF FIGURES 


FIGURE Page 

1 Geometry of vegetation canopy model 5 

2 Layout of the 1974 experimental plots 9 

3 Hemispherical photograph of foliage selections-2 weeks 
before anthesis planted 14 cm apart in equidistant 

spacing pattern 12 

4 Geometry used in derivation of radiative transfer 

equation 16 

5 Geometry used in derivation of general radiative 

transfer equation 18 

6 Geometry of cylindrical scatterer model 24 

7 Sensitivity plot of cylindrical scatterer model with 

different fractional volumes at soil moisture 4.4% ... 29 

8 Sensitivity plot of cylindrical scatterer model with 

different fractional volumes at soil moisture 10.2%. . . 30 

9 Sensitivity plot of cylindrical scatterer model with 

different fractional volumes at soil moisture 25.3%. . . 31 

10 Sensitivity plot of cylindrical scatterer model with 

different e s ' at soil moisture 4.4% 32 

11 Sensitivity plot of cylindrical scatterer model with 

different s s ' at soil moisture 10.2% 33 

12 Sensitivity plot of cylindrical scatterer model with 

different e s ‘ at soil moisture 25.3% 34 

13 Sensitivity plot of cylindrical scatterer model with 

different e s " at soil moisture 4.4% 35 

14 Sensitivity plot of cylindrical scatterer model with 

different e s " at soil moisture 10.2% 36 

15 Sensitivity plot of cylindrical scatterer model with 

different e s " at soil moisture 25.3% 37 

i i i 


LIST OF FIGURES (continued) 


FIGURE Page 

16 Geometry of ellipsoidal scatterer model 39 

17 Orientation of ellipsoid as specified by the Eulerian 

angle of rotation a, $ and y 40 

18 Sensitivity plot of ellipsoidal scatterer model with 10% 

increase of size at soil moisture 13.8% 46 

19 Sensitivity plot of ellipsoidal scatterer model with 10% 

decrease of size at soil moisture 13.8% 47 

20 Sensitivity plot of ellipsoidal scatterer model with 10% 

increase of size at soil moisture 25.3% 48 

21 Sensitivity plot of ellipsoidal scatterer model with 


different fractional volumes at soil moisture 4.4%. ... 49 

22 Sensitivity plot of ellipsoidal scatterer model with 

different fractional volumes at soil moisture 10.2% ... 50 

23 Sensitivity plot of ellipsoidal scatterer model with 

different fractional volumes at soil moisture 25.3% ... 51 

24 Sensitivity plot of ellipsoidal scatterer model with 


different e s ‘ at soil moisture 4.4% 52 

25 Sensitivity plot of ellipsoidal scatterer model with 

different e s ’ at soil moisture 10.2% 53 

26 Sensitivity plot of ellipsoidal scatterer model with 

different e s ‘ at soil moisture 25.3% 54 

27 Sensitivity plot of ellipsoidal scatterer model with 

different e s " at soil moisture 4.4% 55 

28 Sensitivity plot of ellipsoidal scatterer model with 

different e s " at soil moisture 10.2% . 56 

29 Sensitivity plot of ellipsoidal scatterer model with 

different e s “ at soil moisture 25.3% 57 


TV 


LIST OF FIGURES (continued) 


FIGURE Page 

30 Comparison of theoretical and measured soil 

permittivities . 69 

31 Data matching for bare soil at soil moisture 4.3% .... 70 

32 Data matching for bare soil at soil moisture 10.5%. ... 71 

33 Data matching for bare soil at soil moisture 25.9%. ... 72 

34 Data matching for sorghum field at soil moisture 4.4% . . 74 

35 Data matching for sorghum field at soil moisture 8.7% . . 75 

36 Data matching for sorghum field at soil moisture 10.2%. . 76 

37 Data matching for sorghum field at soil moisture 25.3%. . 77 

38 Modified data matching for sorghum field at soil 

moisture 4.4% 79 

39 Modified data matching for sorghum field at soil 

moisture 8.7% 80 

40 Modified data matching for sorghum field at soil 

moisture 10.2% 81 

B1 Geometry used in derivation of scattered fields from a 

cylindrical scatterer 93 


v 


ABSTRACT 


In passive microwave remote sensing of the earth, a theoretical 
model that utilizes the radiative transfer equations has been developed 
to account for the volume scattering effects of the vegetation canopy. 

Vegetation canopies such as alfalfa, sorghum, and corn are simula- 
ted by a layer of ellipsoidal scatterers and cylindrical structures. 

The ellipsoidal scatterers represent the leaves of vegetation and are 
randomly positioned and oriented. The orientation of ellipsoids is 
characterized by a probability density function of Eulerian angles of 
rotation. The cylindrical structures represent the stalks of vegeta- 
tion and their radii are assumed to be much smaller than their lengths. 
The underlying soil is represented by a half-space medium with a 
homogeneous permittivity and uniform temperature profile. The 
radiative transfer equations are solved by a numerical method using a 
Gaussian quadrature formula to compute both the vertical and horizontal 
polarized brightness temperature as a function of observation angle. 

The theory was applied to the i nterpretation of experimental data 
obtained from sorghum covered fields near College Station, Texas. 
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INTRODUCTION 


In recent years, remote sensing techniques utilizing microwave sen- 
sors are being developed for numerous applications. Major emphasis has 
been placed on measuring soil moisture for use in global monitoring of 
renewable resources using satellite based sensor systems. In addition, 
there is considerable interest in using microwave sensors for geologi- 
cal exploration as well as for tracking systems in military ordinance. 
For all of these applications, the effect of vegetation on the micro- 
wave signal is significant and must be understood to utilize the poten- 
tial of microwave remote sensing techniques. The majority of the 
experimental and theoretical work in this area has been in relationship 
to the measurement of soil moisture LI J-L24J . 

Remote sensing, in general, can be classified into two areas: pas- 

sive and active. In passive microwave remote sensing, one measures the 
electromagnetic power intensity naturally emitted by the medium due to 
turmoil agitation. This measurement is specified by a parameter termed 
the radiometric brightness temperature. The brightness temperature of 
a body is defined as the physical temperature that a blackbody must 
have in order to emit the same amount of radiation as the original 
body. In active remote sensing, one provides a source (e.g. radar) of 
electromagnetic radiation and measures the radiation that is scattered, 
usually in the backward direction. 


In remote sensing of a scattering medium such as vegetation, the 
volume scattering has long been recognized as a dominant factor in both 
active and passive cases. Accordingly, the development of a theoreti- 
cal model is essential both in understanding how the physical proper- 
ties of the medium affect the measurements, and in interpreting the 
remote sensing data. 

The volume scattering can be accounted for by modeling the vegeta- 
tion canopy with either a random medium (random medium approach) or a 
homogeneous medium containing discrete scatterers (discrete scatterer 
approach). In the former case, a random medium has an average permit- 
tivity e 1 and its random part is characterized by a correlation func- 
tion with variance s, horizontal correlation length ip and a vertical 
correlation length z. In the discrete scatterer approach, vegetation 
canopy is modeled as a layer of such scattering medium bounded by air 
above and half-space below, and the scatterers are characterized with 
their sizes, permittivities and effective fractional volume imbedded 
inside a homogeneous medium. 

In the theoretical development of the random medium approach for 
passive remote sensing, the volume scattering effect due to medium 
inhomogeneity was first accounted for by Gurvich et a_l_. [25]. Assuming 
the constant temperature profile, they derived the expressions for the 
brightness temperature of a random medium with laminar structure in the 
single scattering approximation. Incorporating nonuniform temperature 
profiles, Tsang and Kong [26] studied the problem with a radiative 
transfer approach. They [27] also derived the emissivity of a half- 
space random medium with a three-dimensional variation. They then 
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solved the radiative transfer equation to obtain the brightness temper- 
atures for a half-space random medium with three-dimensional variation 

[28] . From the Dyson and Bethe-Sal peter equations, they derived 
modified radiative transfer equations for a two-layer random medium 

[29] . With the method of invariant imbedding, Tsang and Kong [30] 
studied the thermal microwave emission from a slab random medium with 
non-uniform scattering, absorption, and temperature profile in the 
vertical direction. Djermakoye and Kong [31] also used the random 
medium approach and the radiative transfer theory to obtain results for 
a N-scattering layer model with a laminar structure. 

For active remote sensing, Stogryn [32] studied scattering by ran- 
dom dielectric constant fluctuations in the low frequency limit using 
the distorted Born approximation. The bistatic scattering cross sec- 
tions for a random medium with lateral and vertical correlations were 
calculated by Tsang and Kong [33]. Fung and Fung [34] obtained the 
bistatic scattering coefficients from a vegetation-like half-space ran- 
dom medium. Zuniga and Kong [35] studied the scattering from a slab of 
random medium using the Born approximation. Then Zuniga et_ £]_. [36] 
extended the result to the second order in albedo to see the depolari- 
zation effect in the backscattering direction. Zuniga et aj_. [37] also 
studied the scattering from a generalized N-layered random medium using 
the Born approximation. 

In developing the discrete approach for passive remote sensing, 
England [38] examined thermal emission darkening caused by randomly 
distributed isotropic point scatterers in a uniform low-loss dielectric 
half-space medium with a radiative transfer approach. He [39] then 
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solved the radiative transfer equation for a scattering layer contain- 
ing Rayleigh particles. Using Mie scattering phase functions the prob- 
lem has been solved for a half-space medium by Tsang and Kong [40]. 
Tsang et al. [41] applied the Mie scattering model to a layer of cloud 
and rainfall. Chang et al. [42] treated the problem of microwave emis- 
sion from a layer of snow and glacier ice at the nadir viewing angle. 

For active remote sensing case, $h i n _et_ _al_. [43] used an iterative 
and a numerical approach to solve the radiative transfer equation with 
a Rayleigh phase function. Using Foldy's and distorted Born approxima- 
tions, Lang [44] calculated the backscattering cross section. Tsang et 
al . [45] applied the radiative transfer theory within the Rayleigh 
approximation to calculate the backscattering cross section of a layer 
of randomly positioned and oriented small ellipsoids. 

The radiative transfer theory has been used extensively in theore- 
tical studies for remote sensing from scattering medium because it is 
simple and includes multiple scattering effects. It starts with the 
radiative transfer equations which govern the propagation of energy 
through the scattering medium. 

In this thesis, the volume scattering effects of vegetation cano- 
pies are studied by introducing a theoretical model that utilizes the 
radiative transfer equation defined for a medium containing discrete 
scatterers as shown in Fig. 1. The model is interpreted in terms of 
microwave emission. Vegetation canopies such as alfalfa, sorghum, and 
corn are simulated by a layer of ellipsoidal scatterers [45] and cylin- 
drical structures. The ellipsoidal scatterers represent the leaves of 
vegetation and are randomly positioned and oriented. The orientation 
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Fig. 1. Geometry of vegetation canopy. 
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of the ellipsoids is characterized by a probability density function of 
Eulerian angles of rotation [46]. The cylindrical structure represents 
the stalks of the vegetation and has a small radius in the horizontal 
direction and a large length in the vertical direction. The underlying 
soil is represented by half-space medium with a homogeneous permittivi- 
ty and a uniform temperature profile. The model computes both the ver- 
tical and horizontal polarized brightness temperature as a function of 
the emission angle. 

The second section of this report describes the experimental mea- 
surements obtained at the Texas A&M University Experimental Research 
Farm near College Station, Texas. The passive microwave and ground 
truth data collected from bare soil and vegetation fields are tabula- 
ted. Also, the information of test vegetation (sorghum) such as the 
shape, size and leaf inclination angle is described. 

In the section entitled Radiative Transfer Theory, the radiative 
transfer equations that describe the electromagnetic emission and scat- 
tering are derived for a homogeneous medium containing randomly dis- 
tributed cylindrical and ellipsoidal scatterers, respectively. For 
each discrete scatterers model, the absorption and scattering coeffi- 
cients, the loss per unit length due to absorption and scattering, 
respectively, as well as scattering functions are derived. Then, the 
radiative transfer equations are solved by numerical technique. 

The next section describes the numerical technique used to solve 
the radiative transfer equations using the Gaussian quadrature formula. 

In the section entitled Data Matching, the theoretical results of 
the radiative transfer equations are matched to the experimental 
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measurements. First, experimental measurements obtained from bare soil 
are used to calculate the permittivity of bare soil. These 
permittivities are used to compute the theoretical results for a 
vegetated soil. These computations are matched to experimental 
measurements of the vegetated fields by adjusting the model input 
parameters such as size, permittivity and volumetric density of each 
scatterer. 
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DESCRIPTION OF MEASUREMENT AND VEGETATION 


Background 

There has been a series of experimental ground-based microwave mea- 
surements at the Texas A&M University Experimental Research Farm near 
College Station, Texas L8J-L12]. An experiment was performed in 1974 
during which simultaneous active and passive microwave measurements 
were acquired over controlled bare and vegetated soil as a function of 
soil moisture and surface roughness. The NASA Johnson Space Center 
L-band (1.4 GHz) and X-band (10.69 GHz) truck mounted radiometer was 
used to obtain the passive microwave data. This thesis deals only with 
passive microwave remote sensing. 

Data Acquisition 

Ground truth and microwave data were obtained from ten plots of 
land in the Texas A&M University Experimental Research Farm in Burleson 
County, Texas [12J. Fig. 2 shows the field layout. Textural analysis 
of soil within the test field revealed that the average soil texture 
was 3% sand, 35% silt and 62% clay. The soil was uniform across the 
test fields and is classified as Miller clay. 

In Fig. 2, nine of ten plots were 15 m by 46 m (sets A, B, C) and 
one was 46 m square (plot D). The nine smaller plots were in sets of 
three with each set having one plot prepared with a smooth surface, one 
with a medium rough surface, and one with a rough surface. Set C was 
bare while sets A and B were uniformly planted in a sorghum hybrid. 
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Fig. 2. Layout of the 1974 Experimental Plots. 




Set A was planted approximately five weeks prior to the field measure- 
ments and set B was planted three weeks before. Plot D was also plant- 
ed in sorghum hybrid, but in row running north and south. 

Passive microwave data were recorded simultaneously at each field 
as a function of soil moisture, incident angles, polarization and fre- 
quency [9]. Passive microwave measurements were acquired at 0°, 20°, 
35° and 50° incident angles for both horizontal and vertical polariza- 
tions at 1.4 GHz and 10.69 GHz. The 1.4 GHz data obtained from bare 
soil (CS field) and vegetation-covered (sorghum hybrid) soil (BS field) 
with smooth surfaces are used in this research and tabulated in Appen- 
dix A. 

Ground truth measurements of soil moisture and temperature profile 
were made simultaneously with the microwave measurements [9], Soil 
moisture samples were obtained at several depths down to 15 cm and are 
tabulated in Appendix A. There are 15 profiles for bare soil (CS 1 - 
CS 15) and 7 profiles for vegetated soil (BS 1 - BS 7). The various 
soil moisture conditions were a result of irrigation and the various 
rainfall events that occurred in that time period. Surface roughness 
and vegetation height were obtained at numerous locations throughout 
the experiment. The smooth surface had an rms surface height of 0.88 
cm, the medium 2.6 cm, and the rough 4.3 cm. The average heights of 
vegetation in fields A and B were 138 cm and 125 cm, respectively, ana 
the volumetric density of vegetation in each of these fields was ap- 
proximately 1 %. 
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Description of Vegetation 


Before the formulation of a theoretical model, the vegetation can- 
opy must be specified. Since sorghum was chosen as the vegetation to 
be modeled, informaton such as shape, size and leaf inclination angle 
for sorghum is needed. 

rlavelka [47] studied the effect of leaf type, plant density and row 
spacing to evaluate the performance of several erect and normal leaved 
grain sorghum selections and hybrids. In his research, two erect leaf 
strains from the World Sorghum Collection, P 407 and SC 170, and a nor- 
mal leaved variety, Caprock, were used. Fig. 3 shows the hemispherical 
photographs taken of the 3 different canopies 2 weeks prior to anthe- 
sis. Examination of the photographs showed that a greater percentage 
of leaf area of the erect leaf selections tended toward a more acute 
inclination, while Caprock exibited most of its leaves in a horizontal 
plane. In Fig. 3, we observe that the shape of sorghum leaf is long 
and thin. Although the length and width of leaf are different with the 
populations of sorghum, those values are within the range of 40-80 cm 
and 4-8 cm, respectively. Therefore, the effective leaf area is within 
the range of 200tt-500ti cm 2 . The stalk diameters of the three sorghums 
were within the range of 10-25 mm. 

Leaf angle measurements were taken by Havelka from the horizontal 
plane at 2 growth stages: 57 days after planting and 103 days after 
planting. Leaf angle distribution measurements made 57 days after 
planting show that leaves of the erect leaf selections became more 
inclined with increased population density. The greatest prepondence 
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of leaf angle measured from the horizontal plane fell into the 45° - 
90° angle classes for P 407 and SC 170. About 70% of the leaf areas of 
P 407 and SC 170 occurred in the 60° - 90° angle classes. Caprock had 
its greatest percentage of leaf area in the 30° - 60° angle classes for 
populations of 129X10 3 , 172X10 3 and 258X10 3 plants/ha. At a population 
of 517xl0 3 plants/ha, 48% of its leaf area occurred in the 60° - 90° 
classes. The leaf area of Caprock was distributed mostly in the 30° - 
75° range of classes. The 0° - 30° angle classes accounted for approx- 
imately 18% of leaf area in each population. At 103 days after plant- 
ing, the leaf inclinations had decreased sharply from those measured 
previously. Leaf angles of the two erect leaf strains were mostly in 
the 30° - 75° angle classes, while the leaf angle of Caprock were dis- 
tributed mostly in 0° - 60° angle classes. 

For the theoretical simulation to be performed later, the length 
and width of leaf were chosen as 64.8 cm and 5 cm which correspond to 
the effective leaf area of 324ir cm 2 . Since th leaf is represented by a 
circular disk, the radius of disk can be assumed to be 18 cm which cor- 
responds to the effective leaf area of 324 rr cm 2 . Furthermore, for the 
convenience of problem, we assume the leaves are uniformly distributed 
along the height of sorghum. 

From the Havelka's study, the leaf angles of erect leaved sorghum 
and Caprock were distributed mostly in 60° - 90° and 30“ - 75’ from the 
horizontal plane, respectively . Thus, a leaf angle of 60° from the 
horizontal plane (30° from the vertical direction) can be chosen for 
the theoretical modeling. 
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RADIATIVE TRANSFER THEORY 


In this section, the radiative transfer equations that govern the 
electromagnetic propagation through the homogeneous medium containing 
discrete scatterers are derived. 

In general, two theories are currently being used to incorporate 
scattering effects: wave theory and the radiative transfer theory. In 

the wave theory one starts out with Maxwell's equations, introduces the 
scattering and absorption characteristics of the medium, and tries to 
solve the problem for the quantities of interest, such as brightness 
temperatures or backscattering cross sections. This is mathematically 
rigorous in a sense that in principle all of the multiple scattering, 
diffraction, and interference effects can be included. However, in 
practice, it is impossible to obtain a formulation which completely 
includes all these effects, and various theories which yield useful 
solutions are all approximations, each being useful over a specific 
range of parameters. 

Radiative transfer theory, on the other hand, starts with the radi- 
ative transfer equations which govern the propagation of energy through 
the scattering medium. It is assumed that there is no correlation be- 
tween electromagnetic vector fields and therefore, the addition of 
intensities is considered rather than the addition of vector fields. 
However, it has an advantage in that it is simple and, more important- 
ly, includes multiple scattering effects. 
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Therefore, in this chapter the radiative transfer equations are 
derived for the homogeneous medium containing sparsely distributed 
cylindrical and ellipsoidal scatterers. Then, the absorption and scat- 
tering coefficients, loss per unit length due to absorption, and scat- 
tering, respecti vely , as well as the scattering function matrix are 
derived for each model of scatterer. In the cylindrical scatterer 
case, the radius of cylinder is assumed to be much smaller than the 
length of the cylinder. In the ellipsoidal scatterer case, the radia- 
tive transfer theory is applied within the Rayleigh approximation. The 
Rayleigh scattering model can be used when the size of the scatterers 
is small compared to the electromagnetic wavelength of interest. To 
solve the radiative transfer equations, the boundary conditions must be 
satisfied. 


Radiative Transfer Equations 

We now derive the radiative transfer equations which are the funda- 
mental equations governing the variation of intensities in the medium 
which absorbs, emits and scatters radiation. 

Consider a small element of cross section da and length ds in the 
medium as shown in Fig. 4. From a phenomenological point of view, 
energy conservation requires that the change in intensity I in a 
distance ds be given by 

dl = - k a Ids - k s Ids + E t ds + k $ Jds (1) 

where k a and k s are the absorption and scattering coefficients. 

The term -k a Ids accounts for extinction due to absorption, the term 
-k s Ids accounts for extinction due to scattering out of the direction 
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of propagation, the term E^ds arises from thermal emission and the 
term k s Jds accounts for the increase in the intensity due to scatter- 
ing from all other directions into the direction of propagation. 

We can generalize the above equation to the case of radiation with 
an arbitrary polarization. Consider a scattering dielectric slab of 
thickness d, permittivity and temperature T, overlying a homogeneous 
half-space of permittivity e 2 and temperature T 2 (Fig. 5). Both the 
scattering layer and the homogeneous half-space are emitting thermal 
radiation. Using the Stokes parameters, the radiative transfer equa- 
tions in matrix form can be written down readily. In the following, we 
take z to be the vertical direction and denote the direction of propag- 
ation by solid angle variable a, where 

fl s ( 9 *<j>) (2) 

Then, 

cos 9 r (n,z) * -? e (a) • r (q,z) + E t 

+ / da' P (a.a 1 ) • I (a'.z)- (3) 

where I (a,z) is a four-component vector which represents the four 
Stokes parameters, and the extinction coefficient, K (a) = Sc. (a) 

6 a 

+ S (a), forms a diagonal matrix. 

£ = Ck v . k h i ^3 « * 4 J (4) 

p(a,a‘ ), the scattering function matrix, represents the scattering 
coupling coefficients from the direction a' = ( © 1 , <j»* ) into the 
direction a = { 9 » 4>) [483 - The coefficients cosa on the left hand side 
of (3) arises from the fact that the propagating direction is inclined 
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Fig. 5. Geometry used in derivation of general 
radiative transfer equation. 
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at angle 8 to the vertical direction z direction, so that ds=sece dz. 
Stokes vector T(fl,z) contains four Stokes parameters which are 


I(a.z) 


l h (a.z) 

U (n,z) 
V (Q.Z) 


(5) 


The first component, I v , denotes the specific intensity for vertical 
polarization and the second component, I&. denotes the specific in- 
tensity for horizontal polarization. The third and fourth components, 
1) and V, represent the correlation between the two polarizations. For 
a time harmonic field given by E = v E v + h where v and h denote 
the two orthogonal polarizations which are perpendicular to the pro- 

A A A A 

pagation direction k, (i.e. v x h - k) , the Stokes parameters are 
defined as 


-r<M 2 > 

(6) 

-k- < M 2 > 

(7) 

Real < E > 

n v h 

(8) 

Imag < E Ef > 
n v h 

(9) 


where n is the wave impedance of the medium and the bracket <> denotes 
the time average. The average is necessary because of the fluctuation 
in the amplitude and phase of the wave. 

The thermal emission coefficient E^in (3) can be derived assuming 

that each point in the medium is in local thermodynamic equilibrium. 

In that case, it is assumed that the circumstances are such that we can 
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define at each point in the medium a local temperature T such that the 
emission coefficient at each point is given by Kirchhoff's law. 

E t = £ a (Q ) B (10) 

where B is the Plank function, 

- hv 3 1 

3 = — r 

X 2 exp[hv/kT]-l 

where K and h are the Boltzmann and Plank constants, v is frequency, c 
is the speed of light in the medium, and T is the temperature of the 
medium. The last two elements of the Plank function £ are zero since 
the emitted thermal radiation is uncoupled. 

In the Rayleigh-Jeans approximation, which is valid in microwave 
frequencies (hv/KT <<1), we obtain from (11) 

8 = CT (12) 

where 


T 

T 



and x is the free space wavelength and ej' is the real part of the 
permittivity of the scattering medium. 

In passive remote sensing case, only two Stokes parameters I v and 
Ih are necessary for the radiative transfer equations. Therefore, 
the radiative transfer equations inside the scattering medium take the 
following form. For 0 < 8 < u 


1 

1 

0 

0 


( 11 ) 
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c°s 0 I v ( e, 4 >,z) = k av ( 0 »4>) CT - k ev ( 9, 4>) I v (e,<M) 


+ J d 0 ' sin 0 ' [(v,v‘) I ( 0‘ , <j>’ ,z) 

0 0 v 

+ ( v,h‘ ) l h (0 , ,* , ,z)] (15) 


cose = k ah ( 8 »4>) CT " k e h(M) V 8 ^* 2 ) 

2 

+ J d0‘ sin 0 1 J d(t>' [(h,v') I ( 0 1 , 4 >' , z) 

0 0 v 

+ (h ,h* ) I h ( e* , 4 >' ,z) J (16) 

where the subscripts v and h represent vertical polarization and 
horizontal polarization, respectively, the extinction coefficient 
k e ( 0 , <{>) = k a ( 0 , 4 >) + k s (e, 4 >) is the loss due to absorption and 
scattering, and the scattering function (v,v‘), (v,h’), (h,v') and 
(h,h') denote the scattering coupling coefficients from the direction 
0', (j, 1 into the direction 9 , <j>. 

Boundary Conditions 


The boundary conditions for the Stokes parameters at planar 
dielectic interface have been derived by Tsang and Kong [33] from the 
continuity of tangential electric magnetic fields. For passive remote 
sensing, the boundary conditions that have to be satisfied by the 
solutions to the equations of (15) and (16) are (Fig. 5), for 0<d<n/2, 
at z = 0. 


I v U-9,Z=0) = 0 

(17) 

I h (n-0,Z=O) = 0 

(18) 
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and at z = -d 


I y ( 6,Z=-d) = f vl2 (3) I v (ir-6,z=-d) + t v21 (9) CT 2 (19) 

l^i (9»2 =_ d) ■ r^i 2 ( 0) I^( 7r-0,z = -d) + t^ 2 i( a) CT 2 (20) 

where 

^v21^ 9 ^ = 1 ■ r v21^ (21) 

t h21 ^ 9 ^ = 1 “ r h21^ 9 ^ ( 22 ) 

and for z, m = 1,2 

r vm< 9) ’ | s tm< 9) | 2 (23) 

W 9 > ■ | R W 9 >I 2 < 24 > 

where R^ and are the TE and TM reflection coefficients, 
respectively, and 


where R^ and are the TE and TM reflection coefficients, 
respectively, and 


V 9 ) ■ 

k 

iz mz 

k. +k . 
zz mz 

(25) 

W 9) ■ 

s m k zz" e z ,( mz 
e m k 2 z +e 2 k mz 

(26) 

with 



k « 

* k 2 C0S8 t 

(27) 

k z 

= w/ u_e. 
H o l 

(28) 

Once the specific intensities 

I v and If, are obtained, the 


brightness temperatures are given 

as 


— 1 

03 

< 

CP 

II 

"Jr 
*— • 

- r vl0 > I v (8. Z=0) 

(29) 


T Bh < 9 ) =4-(l - r hl0 ) I h (e, z=0) (30) 
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Cylindrical scatterer Model 


Formulation 

Consider a collection of sparsely distributed cylinders with 
permittivity e s embedded in region 1 with permittivity e^eo above a 
half space of homogeneous dielectric with permittivity e 2 as shown in 
Fig. 6 . The radiative transfer equations inside region 1 for 0<8<ir 
take the form of equations in (15) and (16). These equations can be 
simplified by assuming that its radius and length are small and long in 
horizontal and vertical directions, respectively. Since the 
configuration and the sources are azimuthal symmetric, the intensities 
are independent of angle <j>. 

Using the above assumptions, we find the specific intensities for 
different 9 are decoupled, and the radiative transfer equations of (15) 
and (16) assume the following forms for 0 < 6 <ir 

cose 3 — I v ( 8 ,z) = CQj ( e) - ©) DI v ( ©,z) + Q 2 (e) l h (e,z) 

+ k av^ ( 9 ) CT (3D 

cos 0 *ai~ r h( 9 » z ) = W2<0) V 9 ’ 2 ) +CQ 3 (e) - kg^( 8)3 l h ( 0 .z) 

+ ^h } (0) CT (32) 


where the absorption coefficients are: 



-4 Na £ __2 

1 — TT Q 

kjcose 



(33) 


k 


(c) 

ah 


( 8 ) 


-4 Na % 
cose 


„a 2 41 


I 

m 


e i _e c 

irrr) 

e l £ s 


(34) 
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Fig. 6. Geometry of cylindrical scatterer model. 
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and the extinction coefficients are: 


k « )(9) ■ k iv } < 9 > + k S^ (e) (35) 

k ih ,(9) ■ k ih > (9) + k sh’ (8) (36) 


with the scattering coefficients, 


k 


(c) 

sv 


( 9 ) 


4 Na £ 
kjcosa 



1 



k 


(c) 

sh 


(9) 


4 Na i 
k x cosa 


[ 2|^ e | 2 * 2 l c > E l 2 ! 


and the scattering functions. 


(37) 

(38) 


where 


with 


Q X ( 9 ) 

= 

4 Na £ 
k L cose 

[|C™p * 2 |C™| 2 ] 

(39) 

Q 2 (9) 

= 

4 Na £ 
kjcosa 

ro 

o 

m 

rs> 

(40) 

0 3 (9) 

3 

4 Na £ 
k^ose 

2 |cT e | 2 

(41) 


-TM 

*-o - 4 

k, 2 

— ( e s*^) 

(42) 


TM iira 2 k lz 2 e $ -l 
C 1 * 4 c s .l 

(43) 


-ME 

'-l 

iira 2 k 1 

r 

k lz V 1 

(44) 


c TE 

L l 

i-n-a 2 k t 
4 

2 V 1 
£ s +1 

(45) 



k lp * 

k x sin a 

(46) 



k lz = 

k L cos a 

(47) 
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In (33) -(41) , N a is the number of scatterers per unit area. These 
coefficients are derived in Appendix B in detail. 


The Solution for Cylindrical Model 

The radiative transfer equations for the cylindrical scatterer 
model can be solved explicitly by obtaining the particular solution and 
homogeneous solution, and by summing them. Then, the solution of 
radiative transfer equations, (31) and (32), subject to the boundary 
conditions in section B are: 
for the upward intensities ( 0<8 <tt/2), 

I v (9,z) = CT + A u e a i z/cose + B u e a 2 z/cos9 (48) 

I. ( 9 ,z) = CT + A e a l z/cos0 + b ^£22. e a 2 z/cose ( 49 ) 

n U L 2 U L 2 

and for downward intensities (tt/2<6<ti), 

I v (9,z) = CT + A d e " a i z/cos0 + B d e“2 z/cos0 (50) 

I. (0,z) = CT + A e a i z/cose + B w ^422. e -a 2 2 /cos 0 (51) 

n u L 2 Q 2 

where the subscript u and d represent the upward direction and downward 

direction, respectively and 

A d = t c 2 "( C l + °2 ) ^ CT/(a 2 -ai) (52) 

= [C 2 -(Ci + ai)] CT/(cq-a 2 ) (53) 

A u = ( f r vl2 " ( C i + a 2^ * r hl2 ( C i + ai ^ A d e ai<i ^ COSQ 

+ ( r vl2 " r hl 2 ^ + a 2 ^ ^d e ° 2 ^ 

/ 

+ [t v i 2 (Ci + a 2 ) - t hl2 C 2 ] (CT 2 - CT)} e a l d/cOS0 / ( a 2 - ai ) (54) 
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8 u ■ 1 (%12 - r hl2> i C l * ».) A d e “ ld/C0S9 
+ [r vl2 (C, + „,) - r hl2 (Cj + « 2 )] B d e a 2 d/cos0 

+ [t v i2 < C i + 01 ) - t hl2 C 2 ] (CT 2 - CT) } e a 2 d/cos9 , (oi . ^ (55J 

with 


°1 = ’[ c l 

+ c 3 

- /" 

(C 1 -C 3 ) 2 + 4C 2 ] / 2 

(56) 

o 

1 

II 

c\J 

a 

+ c 3 

- f 

TCi-CTT 2 ' "+ 4C 2 ] / 2 

(57) 

Cl 

= Qi 

(0) 

- k‘v )(6) 

(58) 

C 2 

= q 2 

(9) 


(59) 

C 3 

- 

(0) 

- k eh )(e) 

(60) 


The brightness temperature for the vertical and horizontal 
polarization can be obtained by applying the boundary conditions at 
z=0. Then, 

T Bu (9) . -J- (A u * B u + CT) (61) 

T Bh (8) =4- [A U (C ! + ai )/C 2 t B U (C 1 + <x 2 )/C 2 + CT ] (62) 

The results 

The formulas derived in this section were used to simulate the 
emissivity for various soil moisture contents at a frequency of 1.4 
GHz. Before matching the measurement data, the sensitivity of the 
emissivity to changes in fractional volume and permittivity of the 
cylindrical scatterers (or stalks) is studied. 
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In Figs. 7-9, the vertical polarized emissivity and the horizontal 
polarized emissivity are plotted as a function of observation angle by 
changing the fractional volume from 0.01% to 0.1% at soil moisture con- 
tents of 4.4% , 10.2% and 25.3% by weight. The permittivity of the 
scatterer is e s 3 (2 + iO.l) e 0 which is an assumed value. The 
radius and length of the scatterers are 1 cm and 125 cm, respectively. 
The permittivities of underlying soil are e 2 3 (3.65 + i 0 . 3 ) e 0 , (8.13 
+ i0.68)e 0 and (27.3 + i2.06)e 0 , respectively and the surface tempera- 
tures are 302. 1°K, 299. 9“K and 292. 1°K, respecti vely. The depth of the 
scattering medium is 125 cm which is the same value as the cylinder 
length. In Figs. 7-9, by changing fractional volume from 0.01% to 
0.1%, emissivity was increased by 20% at the observation angle of 20° 
for each soil moisture content. 

In Figs. 10-12, the sensitivity of the emissivity to the real part 
of scatterer permittivity, g s ‘ , is shown; An increase in e s ’ 
results in a decrease of emissivit' . For the higher soil moisture con- 
tents, emissivities were decreased more than for. the drier, soil mois- 
ture content. 

In Figs. 13-15, the sensitivity of the emissivity to the imaginary 
part of scatterer permittivity, es" , is shown. An increase in e s " 
results in an increase in emissivity. For higher soil moisture con- 
tents, emissivities were increased more. In Figs. 10-15, fractional 
volume was 0.1% and other parameters were the same as those in Figs. 
7-9. 

In the sensitivity analysis, we observe that small changes of 
fractional volume and permittivity of the cylindrical scatterer (or 
stalk) result in relatively large changes in emissivity. 
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different fractional volumes at soil moisture 4.4 %. 
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different fractional volumes at soil moisture 10.2 %. 
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different fractional volumes at soil moisture 25.3 %. 
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Mg. 10. Sensitivity plot of cylindrical scatterer model with 
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Fig. 11. Sensitivity plot of cylindrical scatterer model with 
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My. 12. Sensitivity plot of cylindrical scatterer model with 
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Fig. 13. Sensitivity plot of cylindrical scatterer model with 


Yl. Jill iI'U hL MODEL (change of e 


o 



o 


o 

o 


o 

Q 

o 


o 

o “ 

o 

i7> 

05 

K 

’0 

ill 

T 

•'O 

iij 


o 


o 

o 

o 

3 

O 

Q 

c 

o 

o 

c 


.u:aiss:k2 


36 


fig. 14. Sensitivity plot of cylindrical scatterer model with 
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Kiy. 15. Sensitivity plot of cylindrical scatterer model with 


Ellipsoidal Sc atterer Model 


In this section we apply the radiative transfer approach to 
calculate the brightness temperature of a layer of randomly positioned 
and oriented ellipsoids over a homogeneous dielectric half space. The 
orientation of the ellipsoids is characterized by a probability density 
function of the Eulerian angles of rotation. 

Formulation 

Consider a collection of sparsely distributed ellipsoids with 
permittivity embedded in region I with permittivity e x = e 0 above a 
half space of homogeneous dielectric with permittivity e 2 (Fig. 16). 

The radiative transfer equations inside the region 1 for 0<6<tt take the 
form of equations (15) and (16) with 

k ev (0> + ) ' k lv )(e >*> (63) 

k etl (8,*) * k<jj>( 9lt ) (64) 

where the extinction coefficient is a summation of the absorption and 
scattering coefficients. 

Consider an electromagnetic wave incident from the angular 

direction n' = ( e * » <j> ’ ) on a small ellipsoid with semiaxes a, b, and c 
with the axes pointing in the x^, y^, and z^ directions, respecti vely . 

The orientation of the body of the ellipsoid is specified by the 
Eulerian angles of rotation a, B and y [46] (Fig. 17). 
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Fig. 16. Geometry of ellipsoidal scatterer model. 
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= x (cos y cos a - sin y cos 8 sin a) 

A 

+ y (cos Y Sin a + sin y COS s cos a) 

+ z sin y sin 8 (65) 

A A 

y^ =-x (sin y cos a + cos y cos 3 sin a) 

A 

+ y (-sin y sin a + cos y cos 3 cos a ) 

A 

+ z cos y sin 3 (66) 

A A A A 

Zjj 3 x sin 3 sin a - y sin s cos a + z cos 3 (67) 

To obtain the scattering matrix and extinction coefficients for a 
collection of randomly oriented ellipsoids, we further average over the 
Eulerian angles of rotation a, 3, and y [45]. The derivation of these 
coefficients will be treated in Appendix C. Since the ellipsoidal 
scatterers represent the leaves of vegetation, a good approximation can 
be obtained by assuming the ellipsoid to be a circular disk 
characterized by a 3 b and c << a [45]. In this case, the Eulerian 
angle y is immaterial and the orientation of the ellipsoid can be 
characterized by a and 8. We further assume that the joint probability 
density p(a,s) is uniform in a. Thus 

p(a,S) 3 -2^- 0< a< 2^ (68) 

Then the scattering functions are, for 0 < 9,3' < 17 , 0 < < 2*: 

1^ * 

( v,v ' ) 3 — l*-v 0 f/ ds P(S) {[sin 8 sin 9‘ + cos 8 cos 8’ 

16n 0 

cos ($ -<fr')] 2 | a 9 | 2 + [ C0s2e C0 | 29 ' s — - 8 - (1 + 2 cos 2 (#-♦')) 

. sin 2 8 cos 2 8‘ sin 2 s cos 2 s , 

+ 2 + 
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2 sin 9 ‘ cos 9 ‘ sin 9 cos 9 sin 2 3 cos 2 s cos( $-<♦>* ) 

♦ JL^S-S^e co^j + s1 „2 e sjn 2 8 . co ^ e j . 

j a y ” a ej 2 + 2 [ Sin 9 s-in 9 ’ + C0S 9 C0S 9 ' C0S ( ) ] * 

[ a& 9 cos ^g 81 "* 8 « 8 (♦-♦•) + sin 0 Sin 9* cos 2 a] . 

Re[a* {a^ - a Q )]} ( 69 ) 

(v,h‘) = 1 ~2- v 0 ds p(s) {cos 2 8 sin 2 ( $-$' )| a 0 | 2 + 

[ gifg^t (1 * 2 sin 2 ($-$' ) } + c ° s2 » ] . 

ja^, - a 0 J 2 + cos 2 8 sin 2 g sin 2 ( ) Re [a*(a Y -a 0 )J} ( 70 ) 

(h,v‘) » v 0 ds P(S) {cos 2 9'sin 2 ( '!>“'{>') | a 0 | 2 + 

[ {1 + 2 sin 2 (*-*')) + Sin 2 8| COS 2 a slnVj . 

I a r " a dl 2 + cos2e sin20 sin2 U-<^') Re Ca^(a y - a 0^} (71) 

(h,h‘) » l ^ T v 0 ds P(B) {cos 2 U-$') | a 0 J 2 + 

(i + 2 cos 2 U-*‘))] |a^-a 0 | 2 


+ sin 2 s cos 2 ( ) Re [a*(a Y -a 0 )]} 

(72) 

Lang [44], w§ define 


a _ V 1 

Y l+v d A 3 

(73) 

_ _ v 1 

9 1+v^A 

(74) 

k = w / u 0 e 0 

(75) 


42 


with 


abc 


(e s -l) 


(76) 


IT 

2T 


2 


(77) 


(78) 


In (69)-(72), f is the fractional volume occupied by the disks and, 

f = n Q v 0 (79) 


where n 0 is the number of circular disks per unit volume, and v 0 is the 
volume of one circular disk. 

v 0 ■ yJ-a 2 c (80) 

^(e',*') 3 “5^ v o d S P<8) {| a 9 | 2 + [1/2 cos 2 8'sin 2 0 + 
sin 2 8'cos 2 g] |a^-a g j 2 + 2[sin 2 s'cos 2 g + 

l/2cos 2 8.' si n 2 s] Re [a’U^ -a g )J} (81) 

k sh ) ( 0 ‘ * < * , ‘ ) 3 v 0 f^d8p(8) { j a 9 | 2 + 1/2 sin 2 3 ja y -a 0 | 2 + 

sin 2 s Re U'U^, -a 0 )]} (82) 


and the absorption coefficients are 

kav^ 9 ',*') 3 / ds P( 8 ) [b g + [1/2 cos 2 8‘sin 2 g + 

sin 2 8'cos 2 8) ] (b Y -b e ) ] (83) 

k ih J ( 9* »<4>' ) 3 7 d S P(8) [b e + 1/2 sin 2 s (b Y -b 0 )j (84) 
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where 




b 

7 

= kf 

Imag(a^) 

(85) 



b e 

= kf 

Imag(a Q ) 

(86) 

From (81)- 

■(84), it 

is obvious 

that 

the extinction coefficients 

are 

dependent 

of <j> and 

the following relations hold 




k< e >( 
sv x 

7T — 0 ) 

- k sv )(a > 

(87) 



k< e >( 
K sh 1 

7T— 0 ) 

* k ih’< 9 > 

(88) 



k< e >( 
av v 

1T-0 ) 

- k L e) < 9 ' 

(89) 



k (e), 

ah 1 

1T-0) 

■ k iS , <*) 

(90) 


The Solution and the Results. ' 

The radiative transfer equations in (15) and (16) together with 
the boundary conditions of (17)-(20) can be solved by a numerical 
approach. The numerical approach using a Gaussian quadrature formula 
will be described in the section entitled Numerical Approach. 

In simulating the theoretical model, a probability density function 
p(3) is needed to describe the leaf angle distribution. We shall use 
one of the following four probability density functions for the 
Eulerian angle: 


P 1 ( B ) = 

1 1/ A0 1 1 , 

0 < 0 < AS n 
elsewhere 

(91) 

P 2 (0) = 

r sing/[l-cos( abu)3, 
i 0 

0 <_ 0 <_ AB 1 1 
elsewhere 

(92) 

P 3 (0) = 

rl/2Afri. 
1 0 

tt/2-a&J- _< b tt/2 + ab-L 
el sewhere 

(93) 
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*/2 0 £ ir/2+Afri- /Q 4 * 
elsewhere ' 


Ms) = ( S J" 8/L2 s,n4&ij - 
Probability p x < s) and P 3 (s) were used by Lang [44]. Tsang et_ jal_. [45] 
introduced the probabilities p 2 (s) and p 4 .( s) to take into account the 
smaller angle around 0=0. For a vegetation field, probability pi+(0) 
can be used by assuming that the leaves are inclined more in the ver- 
tical direction. 

Here, like the case of the cylindrical scatterer model, a sensitiv- 
ity analysis of the input parameters were performed. Fig. 18 shows the 
effect of the change of ellipsoidal scatterer (or leaf) size on emis- 
sivity at a soil moisture content of 13.8%. Emissivities are compared 
for two different cases; for the nominal values a, b and c chosen to be 
18 cm, 18 cm and 0.05 cm, respectively and 10% increases of each 
value. The underlying soil was assumed to have a permittivity of e 2 3 
(14.06 + i 0. 95) eg and a surface temperature 291. 2°K. The fractional 
volume of the ellipsoidal scatterers was 1% and leaves were assumed to 
be inclined by 30" from the vertical direction. The depth of the scat- 
tering medium was 125 cm. In Fig. 18, the emissivities were almost the 
same for both cases. For a 10% increase of size, the comparison was 
almost the same as shown in Fig. 19. From the above results, the emis- 
sivity computed by the ellipsoidal scatterer model was not affected by 
a 10% change of scatterer size. The same result was obtained for a 
soil moisture content of 25.3% as shown in Fig. 20. 

Figs. 21-23 show the sensitivity of emissivity to fractional volume 
for different soil moisture contents. 3y changing the fractional 
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Fig. 18. Sensitivity plot of ellipsoidal scatterer model with 
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Fig. 19. Sensitivity plot of ellipsoidal scatterer model with 
10 % decrease of size at soil moisture 13.8 %. 
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Fig. 20. Sensitivity plot of ellipsoidal scatterer model with 
10 % increase of size at soil moisture 25.3 %. 
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different fractional volumes at soil moisture 4.4 %. 
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different fractional volumes at soil moisture 10.2 %. 
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different fractional volumes at soil moisture 25.3 %. 
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Fig. 25. Sensitivity plot of ellipsoidal scatterer model with 
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Fig. 26. Sensitivity plot of ellipsoidal scatterer model with 
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Fig. 27. Sensitivity plot of ellipsoidal scatterer model with 
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Fig. 28. Sensitivity plot of ellipsoidal scatterer model with 
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Fig. 29. Sensitivity plot of ellipsoidal scatterer model with 


volume from 0.1% to 1%, the emissivity was not affected as much as with 
the cylindrical scatterer model. Note that in the case of the cylin- 
drical scatterer model, the fractional volume was changed from 0.01% to 
0.1%. 

The effects of scatterer permittivity at different soil moisture 
contents are shown in Figs. 24-29. An increase in the value of the 
real part from 2 to 3 decreased the emissivity. An increase in the 
imaginary part from 0.1 to 0.2 increased the emissivity as shown in 
Figs. 27-29. For each change of the real part, the change of emissiv- 
ity was almost the same at each soil moisture content. The emissivity 
was more affected by changes in the imaginary part at the higher soil 
moisture content than for the dry soil case. 

Composite Form of Cylindrical and Ellipsoidal Scatterer 

Now we apply the radiative transfer theory to the scattering medium 
containing both of the cylindrical scatterers and ellipsoidal scatter- 
ed with different volume densities. Suppose that the permittivities 
of both the scatterers are the same. Then, the radiative transfer 
equations can be derived by combining the formulations obtained in 
previous sections C and D. 

For 0 < a < it, 

cos 6 V 9 > z > = k iJ ){0) CT - k J t)(9) V 9 ’ z) 

IT 2 TT 

Q 2 ( 9) I „( 8 > z) + / de ’ sin 9' / d^ 1 
0 0 

[ ( v , v * ) I v (0‘,z) + (v,h-) I h (9\z)] (95) 
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COS 8 -fi- I h (8,z) = k^(e) CT - Q 2 (8) I v ( 8 » z } 

It \ IT 2n 

-k> '(e) I h (s,z) + J da* sin 0 ‘ / d<j>' 
n n 0 0 

[(h,v*) I v (8*,z) + (h,h*) I h (8‘,z) ] (96) 

with 

k aJ )(9) 3 k i5 )(0) + k iv )(9) {97) 

k ih ) (0) 3 *$(*) + k£*)(e) (98) 

^ t} (9) • kgy } ( 0 ) - Qj(e) + k^ } ( 8 ) (99) 

k h t)(8) “ k ih )(0) ‘ Q 3(9) + X(«)( 0 ) ( 100 ) 

where subscripts c and e represent the cylindrical model and 
ellipsoidal model, respectively. 

The radiative transfer equations (95) and (96) subject to the 
boundary conditions in section B can be solved by numerical approach, 
and the results are illustrated in Chapter V. 
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NUMERICAL APPROACH 


For the general case, the i ntegro-differential radiative transfer 
equations must be solved by a numerical approach. The radiative 
transfer equations are rewritten letting y = cos 8. Then, following 
Chandrasekhar [48], the integrals in the radiative transfer equations 
are replaced by a Gaussian quadrature [Appendix D] , which is an 
appropriately weighted sum over 2n intervals between the 2n zeros of 
the even-order Legendre polynomials P 2n (v)* We obtain, for i = ± 1, 
±2, ... ±n 


Tfi I »i (2) = K avi CT - K evi 'vf < 2 > * .1 V 2 > 

j=-n J 

* v h j) r hj i 2 >J 


( 101 ) 


»i h 'hi' 2 * ' K ahi CT - K ehl ! vi < 2 >M a jC(Vi.v.) I vj (z) 

j=-n 

+ h i’ h j) Ihj ^ 


( 102 ) 


where 


I gi (z) = I g ( Wi ,z) (103) 

K x 0 i = (104) 

(a^ ) = (a(y.j), g(yj)) (105) 

a and e denote h or v, x denotes a or e, y-j are the zeros of 
P 2n (y), and aj are the Christoffel weighting functions. There are 
2n yj’s and 2n aj’s. They obey the relation 
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(106) 


3 j = a -J 


“i 3 


(107) 


Relations (101) and (102) constitute a system of 4n ordinary 
differential equations with constant coefficients. 

To solve for the homogeneous solutions to (101) and (102), let 

V 2 ’ = I 8i (108) 

Substitute (108) in (101) and (102) with T equal to zero in order to 
determine the 4n eigenvalues of a and the corresponding 4n 
eigenvectors. In matrix notation 


au • I u 3 -K e • I u + F • a • I u + 8 • a • I d (109) 

-ay • ~ -Kg • r^+B*a»I u +F»a»I d (110) 


where I u and I d are two 2n x 1 matrices 



I 

Vl 
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I 

v-1 

• 


• 
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• 
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vn 


v-n 

l u * 

I 

hi 

• 

r d 3 

I 

h-1 

• 


• 

I 

hn 


• 

I 

h-n 


__ - 


_ 


S3 3 

the matrices y, k gJ and a are 2n x 2n diagonal matrices 


(111) 
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and the matrices F and 8 are forward and backward scattering phase 
functions with dimensions of 2n x 2n. 
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II Uu I! CO 


(v ,v ) 

1 1 

(v ,V ) 
n 1 
(h ,v ) 

1 1 

(h ,v ) .. 
n 1 


• (v ,v ) (v ,h ) 

In 11 

. (v ,v ) (v ,h ) 
n n n 1 

• (h ,v ) (h ,h ) 

In 11 

. (h ,v ) (h ,h ) 
n n n 1 


. (v ,h ) 
1 n 

.. (v ,h ) 
n n 
. (h .h ) 
1 n 

(h ,’h ) 
, . n n 


(v ,v )... (v ,v )(v ,h )... (v ,h ) 

1-1 1 -n 1 -1 1 -n 

(v ,v )... (v ,v ){v ,h )... (v ,h ) 

n -1 n -n n -1 n -n 

(h ,v )... (h ,v )(h ,h )... (h ,h ) 

1-1 1 -n 1 -1 1 -n 

(h ,v )... (h ,v )(h ,h ) (h ,h ) 

n -1 n -n n -1 ... n -n 


( 115 ) 


(116) 


From the symmetry relations for scattering phase functions, note that 


=t 

F = F 
=t 

(117) 

8=3 

(118) 


where superscript t denotes the transpose of the matrix. 

We can reduce the number of homogeneous equations from 4n to 2n by 
defining 


I = I + I - (119) 

+ u d 

i = i - I (120) 

- u d 

Adding (109) and (110), 

ou-I - A . I (121) 

and subtracting (109) from (110) 

ay • I + = B • r_ (122) 
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where 


A s + F • a + B • a 
e 


(123) 


fl 3 -K. + F . a - B • a 
e 


(124) 


Combining (121) and (122), we obtain the eigenvalue problem. 


y* 1 . R.y- 1 . A - a 2 T) . I + = 0 


(125) 


where 1 is an identity matrix. Thus is a is an eigenvalue so is -a. 
After the eigenvalue and eigenvector problem is solved, the solution is 


written in the following form 


2n a, z -a.z 

r 3 l { p, r ,e 1 + P , r e 1 } 

+ 2 +1 -2 +2 J 


(126) 


where a^ is the square root of the 2 -th eigenvalue in (125) and 
the corresponding eigenvectors. The constants ? z and P.^ are to be 
determined from the boundary conditions. 

Also, I. is determined from (121). The result is 


2n a.z -a.z 

• 1,1", I., • * * K t l. t • M 


(127) 


where 


r -»= -b : ' 1 - s • h 


(128) 


Finally, from (119) and (120), 


r - h p r e z +p r e M 

U A 1 2 U2 -2 d2 > 

V s l 




( 130 ) 


where 
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( 131 ) 


! di - '-»> < 132 > 

The particular solutions to (101) and (102) are, for i=±l,...,±n, 

- CT (133) 

i^i - CT (134) 

Therefore, the complete solution is 

2n _ a L z _ -a Z 

I u*I.I p t tl ’-* I d 2 e ! + CT < 135 > 

i-1 

2n . a z Z _ -a £ Z 

I d*J 1 ( p t V l tCT < 136 > 

The 4n unknown constants P z and are to be determined from the 
following boundary conditions. 


I d (z=0) = r 10 • I u (z=0) (137) 

I u (z=-d) 3 r 12 • 1^ (z=-d) + t 12 * CT 2 (138) 


where 



r vl01 


0 


0 

r 

vlOn 

r hl01 


r hl0n 


(139) 
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DATA MATCHING 


In this section, the theoretical results from the composite form of 
the cylindrical and ellipsoidal models are used to match the experimen- 
tal data collected from the vegetated field. In data matching, the 
main objective was to find the set of parameters that are consistent 
with the ground truth observations and that also matches the experiment 
measurements. 

Of the several parameters, the most important are the permittivity 
of underlying soil and the fractional volumes. The first step of data 
matching is to find the permittivity of the soil. Following Schmugge 
[49], the permittivity formula is a linear function of soil moisture. 


-9.9 + 1.38 Sm 

for Sm > 11.5 


^ 2.56+ 0.3 Sm 

for Sm< 11.5 

(143) 

-1.44+ 0.185Sm 

for SnP 11.5 


^ 0.06 • Sm 

for Sm< 11.5 

(144) 


where e 2 ' and e 2 " are real and imaginary parts of soil permittivity, 
respectively, and Sm is the soi 1 moisture. 

Soil permittivity also can be calculated from the definition of 
emissivity for bare smooth soil (CS field). 

e( 0) = 1 - r 10 (e) (145) 

The reflectivity r 10 (e) is given in (23) and (24). In this case, the 
imaginary part of the soil permittivity is small compared to the real 
part and is neglected. 
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In Fig. 30, the real parts of soil permittivity computed using mea- 
surements of emissivity and equations (23), (24) and (145) are plotted 
as a function of soil moisture by weight. The permittivity calculated 
from the 1974 experimental radiometric measurements as an observation 
angle 20° are close to the permittivity values computed with equations 
(143) and (144). 

Figs. 31-33 show the theoretical computations of brightness temper- 
ature' using the permittivity based on the bare soil measurements made 
at 20°. These results are well matched to the experimental measure- 
ments from field CS at various soil moistures and the other observation 
angles. In Figs. 31-33, soil moistures are 4.3%, 10.5% and 25.9%, 
respectively, and the corresponding soil permittivities are e 2 = (3.56 
+ i0.39) e 0 , (8.25 + i 0. 98) co and (28.12 + i3.36)e 0 , respectively. 

Another important parameter is the permittivity of the scatterers. 
Since soil moisture content effects to the moisture content in plant, 
it is natural to assume that the permittivity of the scatterers is a 
function of soil moisture. It is shown below that the vegetation per- 
mittivity must be a function- of soil moisture in order to match obser- 
vations to theoretical predictions. 

The permittivity of the scatterers is e 2 - (2 + i 0 . 1 ) e 0 at a soil 
moisture content of 25.3%. This value was obtained with using cylin- 
drical scatterer model only. Emissivities were first computed using 
the composite form (stalks and leaves) of the theoretical model with 
the constant value of scatterer permittivity given above and compared 
to the measurements obtained from the sorghum field (BS field). Figs. 
34-37 show comparisons of the theoretical simulations and the 
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PERMITTIVITY .VS. SOIL MOISTURE FDR BARE SOIL. 



theoretical and measured soil permittivities. 
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experimental measurements for a vegetated field, and for comparison, 
experimental results of a bare field both for soil moisture contents of 
4.4%, 8.7%, 10.2% and 25.3%, respectively. The corresponding soil per- 
mittivities were e 2 s (3.65 + i0.3)e 0 , (7.21 + i0.58)e 0 » (8.13 + 
i0.68)e 0 and (27.3 + i2.Q6)e 0 , respectively and soil surface tempera- 
tures were 301. 2 U K, 301. 3°K, 299. 9 U K and 292. 1 U K, respectively. In all 
cases, the same input parameters are used. For the cylindrical scat- 
tered (stalks), the radius and length were 1 cm and 125 cm, resepec- 
tively. For the ellipsoidal scatterers (leaves), a, b, and c were 18 
cm, 18 cm and 0.05 cm, respectively. The corresponding leaf area was 
324 it cm 2 which corresponds to the length 64.8 cm and width 5 cm of a 
sorghum leaf. The average inclination angle of the leaf distribution 
was chosen as 30° from the vertical direction. This is within the 
range of most populations of sorghum such as P 407, SC 170 and Cap- 
rock. The fractional volume of the cylindrical and ellipsoidal scat- 
terers were 0.1% and 0.9%, respectively. 

Fig. 37 shows the highest soil moisture content case.- The theoret- 
ical results are well matched to the measurement data. The lower soil 
moisture content cases shown in Figs. 34-36 demonstrate that the 
theoretical results are higher than the measured data. The reason for 
this is that the permittivity of the scatterers, e s = (2 + i.Ql)e 0 , 
was chosen at the soil moisture content of 25.3% and not allowed to 
vary as the soil moisture varied. 

Since the change of soil moisture content affects the moisture con- 
tent in sorghum, it is necessary that the permittivity of the scatter- 
ers be changed as the soil moisture changes. From the sensitivity 
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Fig. 34. Data matching for sorghum field at soil moisture 4.4 %. 
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Fig. 35. Data matching for sorghum field at soil moisture 8.7 %. 
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36. Data matching for sorghum field at soil moisture 10.2 %. 
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fig. 37. Data matching for sorghum field at soil moisture 25.3 %. 



analysis, it was shown that an increase of the real part of the permit- 
tivity results in a decraase of emissivity and that an increase of the 
imaginary part of the permittivity results in an increase of emissivi- 
ty. Here we assume that the permittivity of the scatterers is linearly 
dependent upon the soil moisture content by weight. For soil moistures 
of 4.4%, 8.7% and 10.2% which are approximately 1/6, 1/3 and 2/5 of 
25.3%, respectively, the imaginary parts of scatterer permittivity were 
chosen as 1/6, 1/3 and 2/3 of 0.1 which is the value of 25.3% soil 
moisture content case. With this assumption Figs. 38-40 show that the 
theoretical results are well matched to the measurement data. 
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Fig. 30. Modified data matching for sorghum field at soil moisture 
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(j. 40. Modified data matching for sorghum field at soil moisture 10.2 %. 



CONCLUSION 


As an application of passive microwave remote sensing, a theoreti- 
cal model has been developed to study the effect of volume scattering 
from a vegetation canopy. A vegetation canopy was modeled as a homo- 
geneous medium containing a cylindrical scatterer model representing 
the stalks and an ellipsoidal scatterer model representing the leaves. 
Radiative transfer theory with the composite form of each discrete 
scatterer model was used to obtain the brightness temperature from the 
scattering layer. A numerical technique using Gaussian quadrature 
method was utilized to do this. The theoretical results have been used 
to match measured data obtained from a vegetated field. From these 
results, it can be seen that the agreement is reasonable and the model 
can now be used to simulate the effect of vegetation canopies of dif- 
fering height, structure, and volume density. Important observations 
resulting from this work are as follows: 

(1) The soil permittivity calculated from the experimental radiometric 
measurements at an observation angle of 20° are close to Schmugge's 
formula. 

(2) The effect of fractional volume changes in the cylindrical scatter- 
ed (stalks) is much larger than fractional volume changes in the 
ellipsoidal scattered (leaves). 

(3) Emissivity is not affected by a 20% change of leaf area which cor- 
responds to a 10% change of each semi axis of the ellipsoid. 
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(4) An increase of e s ‘, the real part of scatterer permittivity, 
results in a decrease of emissivity. 

(5) An increease of e s '‘ , the imaginary part of scatterer permittiv- 
ity, results in an increase of emissivity. 

(6) The permittivity of the scatterer is assumed to be linearly depen- 
dent upon the soil moisture, which means that for dry soil, mois- 
ture content in vegetation also reduces and then, reduces the per- 
mittivity of vegetation. Therefore, it is necessary to develop the 
permittivity of vegetation as a function of soil moisture. 

The task of developing a theoretical model of vegetation is by no 
means complete. We have assumed that the underlying soil is a homoge- 
neous half-space. This can immediately be extended to consider the 
case of a multi-layered medium below the vegetation. 
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APPENDIX A 


Table 1. Ground Truth Data for Bare Soil. 


Data 

Set 

Soil Moisture by 
Weight (%) 

Soil Temperature 

(°K) 

CS 

1 

3.1 

313.4 

CS 

2 

4.3 

288.4 

CS 

3 

27.2 

300.5 

CS 

4 

21.8 

304.5 

CS 

5 

19.3 

296.0 

CS 

6 

20.8 

300.9 

CS 

7 

25.9 

304.9 

CS 

8 

24.5 

298.9 

CS 

9 

20.2 

305.8 

CS 

10 

12.8 

299.9 

CS 

11 

15.1 

302.1 

CS 

12 

12.6 

299.5 

CS 

13 

9.5 

302.1 

CS 

14 

10.5 

302.3 

CS 

15 

l'J.4 

300.2 


Table 2. Ground Truth Data for Sorghum Field. 


Data 

Set 

Soil Moisture by 
Weight (%) 

Soil Temperature 

TO 

8S 

1 

i 

25.3 

i 

292.1 

BS 

2 

13.8 

291.2 

BS 

3 

11.7 

303.9 

BS 

4 

10.2 

299.9 

BS 

5 

8.7 

301.3 

8S 

6 

4.8 

310.1 

BS 

7 

4.4 

302.1 
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able 3. Brightness Temperature as a function of 
Observation Angle, a, for Bare Soil. 
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APPENDIX B 


Derivation of scattering functions, scattering coefficients and 
extinction coefficients for a homogeneous medium containing cylindrical 
scatterers. 

Consider a plane wave with the electric field 


E'(r) = (v'E v + h* E h 1 ) e 1 ^‘* r (Bl) 

incident on a long cylinder of radius a, length i and permittivity 

A A 

e s , as shown in Fig. Bl. The unit vectors h‘ , and v‘ denotes the 
direction of the electric fields for horizontal and vertical 
polarizations. The unit vectors h' , v' and k' form an orthogonal 
coordinate system. 


h' * (z * x k' )/| z* x k * J 

AAA 

v' * h* X k' 

In spherical coordinates 

k‘ 3 x sin e‘ cos + y sin 0 ' sin $'+ z cos -a 1 
h' = -x sin <p' + y cos 

A A A A 

v* =» x cos 8* cos + y cos e* sin - z sin 0 ' 
The scattered wave takes the form 

e = ; e v * h ‘ £, 

where 

A A A A 

k = x sin 0 cos <p + y sin 0 sin $ + z cos 0 

A A A 

h = -x sin $ + y cos $ 

A A A A 

v = x cos 8 cos $ + y cos 0 sin $ - z sin 0 


(B2) 

(B3) 

(B4) 
(B5 ) 
(B6) 

( B7 ) 

( 88 ) 

( B9 ) 
(810) 
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The amplitudes of scattered field are related to those of incident 


field by 



i k , ' r 
e 

f 

f 


E'~ 

vv 

vh 


V 

r 

f 

f 


E* 


hv 

hh 


h 


U -J 


— — 


(Bll) 


By the solution of homogeneous Helmholtz equation in cylindrical 
coordinates, the longitudinal components of electric and magnetic 
fields internal to the surface of cylinder take the form of 


E* nt .?A J (k pje 1 ”'*-*') e 1klzZ 
z L m m Sp 


(B12) 


(ip* 


Hl nt =l a j (k p ) e 1 " 12 " 

z m m So 
rrp — 


( B13) 


where J m ( k Sp P ) denote the 8essel function of first kind of order m 
and. 


k s P 

" / k s‘ " k lz' 

( B14) 

k lz 

= k x cos 8 

( B15 ) 

k s 

= u / Uo e s 

(B16) 

k l 

= u / uo £l 

(817) 

Similarly, the longitudinal components of the scattered 

fields take the 


form of 

h ” I C m H m I) { k lo 0> e''” 1 *'*’ 1 e’ klzZ (518) 

IIP— P 

H 2 = f ° m H P > (k p) e’” 1 *-*' 1 e 1kuZ (819) 

(IP— P 
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where H^|k, p) denote the Hankel function of first kind of m and, 
m ip 

k lp = 9 (B20) 

The scattered wave can be solved from Maxwell's equation by 
applying Huygen's principle [50] and by matching the boundary 
conditions. The Huygen's principle is such that with a knowledge of 
the tangential fields over a surface due to original sources, the 
fields everywhere external to the surface can be calculated through the 
use of equivalent new sources. Under the assumption that the 
tangential electric and magnetic fields are known on the surface s' 
enclosing a radiation source, the scattered electric field E(r) outside 
s' takes the form: 

E(r) * $ ds' { iuu G(r,r' ) • [n x H(r‘ )] + v x G(r,r’ ) .[n x E(r' )](B21) 
s' 

where G(r,r') is the dyadic Green's function [51] and, 

§(?,?' ) . (f t 77) (B22) 

With the far field approximation, the amplitudes of scattered field are 
related to those of incident field by 


with 



A = l C™ e 1ffl( +-* ,) • F( 8, 9* ) 


(B24) 


B « ? C* ME e 1m( *“*' } • F( 9,9' ) 


(825) 
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C = l C TE e 1m( *“*' } • F( 9 , 8‘ ) 


C ™ ■ -(V™ - J ,»<V> H i U < k l,*> (WM*!'*. 

- MW - q 2 m "i" < k lp a > C J n,<V>) 2 l' X n 


q m [J m (k c a)] 2 
m m sp 


'm irk. 4 a 

lp 


(B27) 


(B28) 


(B29) 


sin [k 1 |> (cos 0 ' -cos 8) ] 
F( 0,0 ) 3 i <1 ^ ( cos 9 1 -cos 0 ) 


( B30) 


X s * P m N - q 2 [J_(k a)] 2 Qij ' 1 ' 1 (k, a)] 2 

m m m m *- m k s p 7 J ^ m lp 


P m = k i {lc7“ H m 1)(k lp a) J m' (k s P a) * TT H m ) ' (k lp a) J m (k s P a) 1 
Sp lp 

= l^lz r 1 . 1 | 

a I k 2 " k. 2 * 

Sp ip 

M 3 k, (pi— J (k a) J ' (k a) - pi J ' (k. a) J (k a) 1 
m 1 1 k m v p 7 m v s p ' <7 m ' 1 p 7 m's p ,J 

sp ip 


N " = n ( ^; H " 1)(kl » a) J ”' (ks ° a) •^” 1),<k i» a) J ™ ( v>> 


(B31) 


v m w> V‘V> - vV<V> J m (k s 0 a n < B32 > 

i S Q lp 


(B33) 


(B34J 


( B35) 


(836) 


We now derive the scattering functions. Define the radial Poynting 
vector as follows: 


96 


s v=-Hm : 




(837) 

(B38) 


The Stokes parameters are related to the radial Poynting vector by 


! V 3 COS 8 aV r s v 


Iu 


- 1 r Z 3 

cos 0 AV ~ ^h 


(B39) 

(840) 


where 9 is the angle between the direction of propagation and the 
outward normal to the elementary area aA. We then use (B23) to write 
down the components of radial Poynting vector for the scattered field 
in terms of incident Stokes parameters. 


S 

v 





S 

h 




(B41) 


(842) 


Let us now consider a plane wave with the electric field vector 
given by (Bl) incident on a volume aV of homogeneous medium containing 
n cylindrical scatterers with radius a, length i and permittivity e s . 
Relating Stokes parameters of the scattered field to those of incident 
field by making use of ( B39 ) -( 842) , we obtain 


I 



ri A i 2 

U 

CD 


I ‘ 

V 

n 

l 2 

1 1 

i i 


V 

I 

h 

3 

AV 

* 2 

j 8 r 

l c i 2 _ 


I 1 
h 


( 843 ) 
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where we made use of the fact that aV = Aas. From the above expression 
we can easily identify scattering functions 


(v»V* ) 

ir 2 cose‘ 1 
2 

A l 

| 2 

(B44) 

(v,h‘) 

(h,V) 

N a ir 2 cose' 1 
= ( v,h‘ ) 

Bj 

1 2 

(845) 

( 846) 

(h.h*) 

- fi ll 

rt a ir 2 COS0‘ 1 

l C ! 

1 2 

(847) 


where N a in the above equation is the number of scatterers per unit 
volume. 

Once we have obtained the scattering function matrix, we can 
calculate the scattering coefficient k s , which is the functional loss 
of power per unit length due to scattering. 

* sv 3 / do' C(v,V) +(M*)3 (B48) 

k sh = / da* C(v,h' ) +(h,h‘ )] (B48) 

The extinction coefficient which is the fractional loss of power 
per unit length due to absorption and scattering per scatterer, can be 
found by using the optical theorem. The rate of extinction is by 
definition of optical theorem, 

W = - % Imag [E‘ E* k' . k ] (B50) 


where E* is the complex conjugate of the scattered field E. Since the 
extinction power is given by W«aA»N and incident power is aA|E'J 2 cos 0 / 


k (a*) = 
ev' ' 


“V 


ITIcoieLl R e< C ™> 
* m=- 


( 851 ) 


u — TC 

k . (8*) = - t- 2 r V R (C fc ) 

eh' ' k,cose‘ L e v m ' 

1 FTP—® 


( 852 ) 


The radiative transfer equation for cylindrical scatterers can be 
simplified considerably by applying the small radii approximations 


(i.e. k L a<< 1) into (B24)-(B26). 

A(fl.a') = F( e , 6 ' ) CC™ + 2 cj M cos (*-*•)] 
B(a.n’) - Fte.a 1 ) 2 i cf sin (*-*•) 

3 F( 9 , a' ) 2 Cp cos (*-$')] 
with C™, cj M f C? E . and Cp expressed in (42)-(45). 


( B53) 
( B54) 


( B55) 


We can also simplify (B30) with the long cylinder approximation as 
follow: 


lira / d 9 1 sin 3' F 2 (e,9') * ~- 
0 1 

And with the same approximation in (B56) 


(B56) 


/ de * sin 9* F 2 ( 8,8* ) f( 8* ) = rr f(9) 
n i* 


(B56) 


where f(d) is any smoothly varying function of 9. Using these 
approximations, we can find that the specific intensities for different 
9 are decoupled and the radiative transfer equation assume the form of 
(31) and (32) with the coefficients given in (33) -( 38 ) . 
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APPENDIX C 


Derivation of scattering phase function, scattering coefficients and 
absorption coefficients for ellipsoidal scatterers. 

Consider an incident plane wave 

£' = [v* E y * + h* E h ‘ ] e lE i'* r (Cl) 

impinging on an ellipsoidal scatterer of permittivity e s = e s ' + i 
es" embedded in a medium of permittivity ei. The amplitude of the 
scattered field are related to those of the incident field by 


R (a)”| 

ik, r 

f (a,a‘) f . (a.a 1 ) 


V (a*) 

V 

e 

w vh 


V 

E (si) 

r 

f. (a.a*) f (n,n‘) 


E* ( a* ) 

h 


hv hh 


h 


(C2) 


If the axes of the ellipsoid are oriented in the x^, y fa and z b direc- 
tions, then the internal induced dipole moment p 1nt inside the 
ellipsoid is [52] 


-i nt 
P 


(e -e ) 
S 1 


l " Vb 

1 + V* 


E y 

yb b 

+ 

l+u A2 
d 


E z ~ 
zb b 


Uu / 3 


(C3) 


where v 0 is the volume of the ellipsoid, and 


v d " ~T v V e l } 


(C4) 



(C5 ) 



(C5) 
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R ■ [ (s + a 2 )(s + b 2 )(s + c 2 ) ] 1/2 
In (C3) 

E xb 3 < ; *-V E v 

A A 

E yb = (v ' *^b> E v 
E zb = E v 


ds 

(s+c*)k; 


(C7) 



(C8J 

A A 

+ (h‘.x b ) 

E h' 

(C9) 

+ (h‘ *y b ) 

E h‘ 

(CIO) 

A A 

+ (h* -Z b ) 

V 

(Cll) 


The equivalent induced current is 3 n{ * = -iwp 1nt 5(r). This is an 
equivalent Herzian dipole, the radiated fields of which can be 
calculated easily, is in the form of (C2). The coefficients 
Wn.fl’). f V h(fl.a‘), and f hh(fl*fl') are 9 iven by 


f vv (n.fl') 



r .tvx^HXb-V) 
^ 1 + U (j A l 


+ 


A A A A 

(v-y b )(y b »v' ) 
1 + u d A 2 • 


+ 


(v.z^Ub-v 1 ) n 
1 + u d A 3 J 


( C12) 


2 ( v-x. ) (x. -h* ) (v.y )(y -h') 

f vh (n * fl,) 3 17 V ° {e s' s l j C-T+ u'Ai + ~ 1 + ” d A 2 — 


(v-* b )(z b -h- ) n 
1 + u d A 3 J 


( C13) 
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u 2 n ( h -x h )(x h .V) (h-y.) ty h -V ) 

f hv (n,n') * -£ v o (e s -ei) C 1 Vu"A'i + 1 ~ o d A; — 


^ (h.z b )(z b .V) ^ 

+ 1 + Vs ** 


( C14) 


f hh (n,Q') ■ 


(h.x b )(x b .h‘ ) (h.y b )(y b *h‘ ) 


tS v » ‘v««> [ 1 +V° * 


1 + ^ d A 2 


(h-2 b ) (Z b -h* ) 

1 + u d A 3 


■] 


(C15) 


In (C12)-(C15), v‘ and h ' are, respectively, the vertical polarization 

A A 

vector and horizontal polarization vector for incident wave and v and h 
are, respectively, the vertical polarization vector and horizontal 
polarization vector for scattered wave as ^hown in Appendix B. 

Then, with the radial Poynting vectors defined in (B37) and (B38), 
and the relation between the Stokes parameters and the radial Poynting 
vectors as shown in (B41) and (B42) , we can ojtam the scattering 


matrix for passive remote sensing. 

( v, v 1 ) = n 0 < jf yv (n,n')| 2 > (C16) 

(v,h‘ ) = n 0 < jf yh U,fl’)| 2 > (C17) 

(h,v‘ ) = n 0 < |f hv (a,n‘ )| 2 > (C18) 

(h ,h‘ ) = n 0 < |f hh (a,a')| 2 > (C19) 


The angular bracket in ( C16 ) -( C19 ) stands for ensemble average over the 
Eulerian angles, and n 0 stands for the number of ellipsoids per unit 
volume. 
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The scattering coefficient which are the fractional loss of power 


per unit length due to scattering are 

K sv («*) = / da C(v,V) + (h ,v‘ )] (C20) 

k sh (a 1 ) = / da [(h.v 1 ) + (h.h 1 )] (C21) 


The absorption coefficients is n 0 times the fractional loss of 
power due to absorption per scatterer and defined as 


k (a‘) 3 n 0 


/ dv e-'ajjE 1 ^! 2 




2n 


(C22) 


With this definition, the absorption coefficients for ellipsoids are 

„ / (x b .v*) 2 (y b *v') 2 

k av (il '> * ■'If"- v o V " <[ 


♦ ‘ v, ' )2 n 


(C23) 


k ah (n,) 3 n o v o e 5 " w ([ 


(x b *h’) 2 (y b *h’) 2 

1+y dM + | 1+u d A z' 


+ 



(C23) 


where the angular brackets represent ensemble average over the Eulerian 
angle of rotation and can be obtained by multiplying with the probabil- 
ity density function p(o,3.y) and integrating over a between 0 and 2it, 

3 between 0 and 2 tt and y between 0 and Z-a. 
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APPENDIX D 


Gaussian Quadrature Formula 


Consider an integral 


I = / dx f(x) W(x) 
a 


(Dl) 


to be evaluated approximately by using m values of f(x) in the interval 
(a ,b) : 


/ dx f(x) W(x) * l a. f(xj (D2) 

a j=I J J 

where W(x) is a weighting function. The problem is to choose xj and 
determine aj in such manner that the integral is best approximated by 
the least number of terms (smallest m). 

Suppose the m values of f(x) are taken at the positions 

(D3) 


Using the Lagrange interpolation formula we can construct a polynomial 
0(x) of a degree less than or equal to m-1, which will take the same 
values at the points [xj] as f(x) does. Thus 

(»-»ji x p‘(» J ) (D4) 

where 


m 

F(x) = n (x-x i ) (05) 

j-1 J 

is a polynomial of degree m whose zeros are at [xj], and 

F'(Xj) - [-lfF(x) ] .j (Xj-x,) ( 06 ) 
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We define the integral 


b 

lx - / dx W(x) 0(x) (D7) 

a 

Substituting (04) into (07) 


where 



(D8) 


. - 1 b r dx F(x)W(x) 

'i l 


(09) 


Clearly, if f(x) is a polynomial of degree less than or equal to m-1, 
then 


and, therefore 


0(x) = f(x) 


( 010 ) 


l l = I (Dll) 

We note that (D10) and (Dll) hold irrespective of the choices of 


[Xj]« 

We shall now show that with a judicious choice of [xj] 


II 3 1 

(012) 

when f(x) is a polynomial of degree 2m-l or less. 


Let f(x) be a polynomial of degree 0n-l. Since, 

for j*l, . . . ,m 

f(Xj) = 0(Xj) 

(013) 

F <»j) - 0 

(014) 

we can write 


m-1 

f(x) = 0(x) + F(x) l C x l 

0 £ 

(D15) 
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where C a are constants. Since F(x) is a polynomial of degree m, the 
series in the second term in (D 15) is a polynomial of degree 2m-l. By 
proper choice of C^, any function which can be described by a 
polynomial of degree 2m-l can be represented by (D15). Substituting in 
(D15) into (01) 


m m-1 b 

I • l a. f( x . ) + l C / dx F(x) W(x) x* (D16) 

j=l J J t=0 1 a 

Now choosing [xj], and thus F(x), such that 
b 

/ dx F(x) W(x) x l - 0 l 3 0,...,m-l (D17) 

a 

then the second term on the right side of (D16) vanishes and 

I-I x (D18) 

If the polynomial F(x) can be determined from (D 17) , then xj are the 
zeros of F(x) and aj can then be determined from (09). 

For the integrals in the radiative transfer equations W(x) » 1, 
a = -1, b = +1. It is also recognized that the Legendre polynomial 
P m (x) is orthogonal to all x^, and therefore choose 

F(x) = P m (x) (D19) 

and x should be the zeros of the Legendre polynomial P m ( x ) . 

Thus 


where 


m 


/ f(x) dx = l a. f(xj 
-i j=l J J 


1 \ P m< x > . 

J = .{ dx 


(D20) 


( 021 ) 
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This is called Gaussian quatrature formula. The numerical technique of 
using Gaussian quadrature formula to replace the integral in the 
radiative transfer equations and solving the resulting system of first 
order differential equations is called the method of Gaussian 
quadrature. 

For our problem, it is convenient to choose an even order, m=2N, 
Legendre polynomial . Using y instead of x. 



and for j * 1,...,N 



(D23) 

1 } P 2N (w) , 

(D24) 

■ a -j 

(025) 


(026) 


The values of aj and yj for different values of N are tabulated and 
readily available [53]. 
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